A continuous source of translationally cold dipolar molecules 
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The Stark interaction of polar molecules with an inhomogeneous electric field is exploited to select 
slow molecules from a room-temperature reservoir and guide them into an ultrahigh vacuum cham- 
ber. A linear electrostatic quadrupole with a curved section selects molecules with small transverse 
and longitudinal velocities. The source is tested with formaldehyde (H2CO) and deuterated ammo- 
nia (ND3). With H2CO a continuous flux is measured of « 10 9 /s and a longitudinal temperature 
of a few K. The data are compared with the result of a Monte Carlo simulation. 

PACS numbers: 33.80.Ps, 33.55.Be, 39.10.+J 



The past years have seen an explosion of activity in the 
field of cold atomic gases jjj. It is interesting and desir- 
able to extend these investigations to molecules, which 
have a complex internal structure and can as a con- 
sequence possess a permanent electric dipole moment. 
Trapping cold polar molecules will lead to new physics 
due to the long range and anisotropy of the dipole-dipole 
interaction. Slow molecules for precision measurements 
or interferometry are further motivations behind the on- 
going efforts. However, the complexity and density of 
energy levels in the rotational and vibrational manifolds 
largely precludes the effective use of laser cooling tech- 
niques H. Therefore, a number of different approaches 
has been considered for cooling and trapping molecules. 
Buffer-gas cooling in a cryogenic environment is one pos- 
sibility, but requires a rather complex setup || . Another 
method is photoassociation, but this is limited to simple 
molecules with laser-cooled precursor atoms B . A novel 
technique uses deceleration by the Stark effect, where 
packages of polar molecules are decelerated with time- 
varying electric fields || [| [?J. Other, mostly mechani- 
cal methods have also been proposed but remain to be 
demonstrated (|, ||] . 

It is, however, not necessary to produce slow molecules, 
as they are present in any thermal gas, even at room tem- 
perature. Slow molecules only need to be filtered out. For 
this reason, already in the 1950s it was attempted to se- 
lect the slowest atoms from a hot beam using gravity [jl0| . 
These attempts failed, mostly because the slow particles 
were kicked away by the fast ones. Much later, it was 
demonstrated that slow lithium atoms can be efficiently 
guided out of a hot beam with strong permanent mag- 
nets, providing a robust and cheap source of slow atoms 
jnj, e.g. for Bose-Einstein condensation experiments. In 
the same spirit, an efficient and simple filtering technique 
could play an important role towards the production of 
a cold molecular gas. 

In this Letter we describe an experiment in which the 
Stark interaction of polar molecules with an inhomoge- 
neous, electrostatic field is exploited to efficiently select 
and guide slow molecules out of a room-temperature 
reservoir into ultrahigh vacuum. Whether a dipolar 
molecule is weak-field seeking and trapped by an elec- 
tric field minimum, or strong-field seeking and expelled, 



depends on whether the average orientation of the rotat- 
ing molecular dipole is antiparallel or parallel to the local 
electric field, respectively. Weak-field seeking molecules 
are trapped if their kinetic energy is less than the Stark- 
potential barrier for them. Thus a static quadrupolar po- 
tential forms a two-dimensional trap for molecules with 
small transverse velocity components with respect to the 
quadrupole axis. In a linear guide, transverse and longi- 
tudinal components of the velocity are completely decou- 
pled. The longitudinal velocity of the molecules is lim- 
ited by guiding them around a bend in the quadrupole, 
downstream from the reservoir. The centripetal force 
due to the electric field gradient guides only the slowest 
molecules around the bend. The fast molecules escape 
the guide and are pumped away. 

Important parameters for characterizing such a contin- 
uous source of slow molecules are the velocity distribution 
and the flux of the molecules. For a given molecular state, 
we have a maximum transverse (ftmax) and longitudinal 
(uimax) velocity that will be guided. Under the assump- 
tions that the Stark shift of the molecule is linear with 
electric field and that the longitudinal (vi) and transverse 
(i>t) velocities are much smaller than the mean thermal 
velocity inside the reservoir, the flux $ in the guide is 
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where E is the depth-determining electric field strength. 
Eq. (Q) is valid for every molecule with a linear Stark 
state, and therefore also for an ensemble of molecules 
with different, but linear, Stark shifts. The derivation of 
Eq. |l]) utilizes a unique property of the guide, namely, 
that the longitudinal velocity distribution has a linear 
dependence on velocity for small velocities: cx 
vi exp[— vf /a 2 }, where a is the most probable velocity for 
molecules in the source, in contrast to molecular beams, 
which have a vf exp[— vf /a 2 ] dependence. The reason 
is that the guide selects on kinetic energy, not on angle 
as in a typical molecular beam, which is collimated by 
apertures. As a consequence, molecules that are slow in 
axial and radial direction will still be guided, whereas 
they would be lost in a molecular beam. 

Most of our experiments are performed with formalde- 
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FIG. 1: (a) Schematic of the experimental setup, and (b) 
a contour plot of the quadrupolar electric field in a plane 
perpendicular to the electrode rods. 



hyde (H2CO), which has a relatively large dipole moment 
of 2.34Debye and is one of the best studied 4-atomic 
asymmetric-top molecules. At room temperature (300 K) 
almost all the H2CO molecules will be in the vibrational 
ground state. A large number of rotational states consis- 
tent with the Boltzmann distribution at 300 K are pop- 
ulated. The experiment itself (Fig. |l](a)) consists of the 
bent quadrupolar guide starting at a 0.5 mm diameter ce- 
ramic nozzle of 11 mm length and ending in a UHV detec- 
tion chamber, passing through a vacuum chamber where 
the filtering is performed. The guide has a length of 18 cm 
and is made of 1 mm diameter stainless steel rods, with a 
1 mm gap between neighboring rods. The rods are built 
around the ceramic nozzle. Typical operation pressures 
in the nozzle, which injects the molecules directly into 
the quadrupole, are below 0.05 mbar in order to main- 
tain molecular flow conditions. Most of the molecules are 
not guided and escape into the vacuum chamber, where 
a typical operational pressure of a few times 10 -7 mbar 
is maintained by a 500 l/s turbo-molecular pump. The 
90° bent section of the guide is about 60 mm downstream 
from the nozzle and its radius of curvature is 13.5 mm. 
After the bend, the guide passes through a 5 cm long and 
4 mm narrow tube for differential pumping before enter- 
ing the detection chamber, where a 56 l/s turbo pump 
maintains a pressure of » 2 x 10~ 9 mbar. The guide 
abruptly ends 25 mm before the ionization volume of a 
quadrupole mass spectrometer (QMS) [Pfeiffer vacuum, 
Prisma QMA 200] . Apart from measuring the direct flux, 
the QMS monitors the residual gas in the detection cham- 
ber. Standard pressure gauges are used to monitor the 
pressure in the UHV chamber. 

To interpret experimental results, a Monte Carlo sim- 
ulation was performed. The electric field in a 90°-bent 
quadrupole with short straight sections on either side was 
calculated numerically in 3 dimensions. The input distri- 
bution of the molecules injected into the guide in the sim- 
ulation is random and consistent with the thermal speed 
distribution. Its angle distribution is that of a 3 mm deep 
field- free nozzle [ |12| . This is a good approximation to 
the experimental situation, where all mm long nozzle in 
the presence of an electric field results in the transverse 
acceleration of molecules, giving a wall-to-wall collision- 



free path length of approximately 3 mm for the guidable 
molecules. For each molecular species, the Stark shifts 
of the relevant states were calculated by direct numeri- 
cal diagonalization of the Stark Hamiltonian |l3] . In the 
relevant electric field range (0-100 kV/cm), the majority 
of thermally populated rotational states of H2CO and 
ND3 show linear Stark shifts. These were gathered in 
a small number (20) of groups with similar Stark shifts. 
The trajectories under the influence of the electrostatic 
forces were calculated by a Runge-Kutta method. The 
simulation confirmed the expected quadratic dependence 
of the flux on the electric field, see Fig.|](d). The sim- 
ulation also yields the overall transmission of the guide, 
the velocity distribution inside the guide, the output di- 
vergence and velocity distribution of the beam. Finally, 
the simulation shows that the guide enriches the gas with 
states having a large Stark shift. 
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FIG. 2: (a) The high voltage applied to the guide electrodes 
and (b) the molecular flux as a function of time. Note the 
delay and the relatively slow rise of the signal, corresponding 
to a gradual build-up of molecular flux, compared to the sud- 
den fall due to the loss of molecules when the field is switched 
off. The line is a 20-point running average, (c) The flux as a 
function of applied electrode voltage. The symbols denote the 
measured height of the step function in (b). (d) The number 
of trajectories intersecting with the detector in a simulation 
where 10 6 particles were injected into the guide, correspond- 
ing to a 1% fraction of the full input distribution. The lines 
in (c) and (d) are quadratic fits. 



For detecting H2CO molecules the quadrupole mass 
filter is set at mass 29, the strongest peak in the H2CO 
mass spectrum. The influence of other gases at this mass 
is negligible. The channeltron-amplified QMS ion cur- 
rent can be tapped directly for transient measurements. 
Switching on and off the quadrupolar field resulted in a 
modulated QMS signal, see Fig. |(a) and (b). It was 
checked that the electric disturbance from switching the 
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high voltage did not significantly influence the signal 
apart from a short spike of sub micro-second duration. 
At high background pressures of polar molecules in the 
filter chamber and without injection of polar molecules 
into the guide, the ion signal decreases slightly when the 
HV is turned on. This effect is attributed to a decrease of 
conductance of the differential pumping section for polar 
molecules in the presence of electric fields. The effect is 
absent for non-polar molecules. As a further test that the 
signal represents the direct flux of the guide, a mechani- 
cal shutter was installed that can prevent the direct flux 
from reaching the QMS. It is designed so that it only 
prevents the guided molecules from entering the QMS. 
With the shutter blocking the direct flux into the QMS, 
the modulation of the QMS signal was a factor of ks 5 
smaller. With the above tests we conclude that the ob- 
served changes in the QMS signal are due to real changes 
in the density of the guided gas. 
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FIG. 3: The longitudinal velocity distribution (data points 
with statistical error bars) derived from data obtained 
with an electrode voltage of ±5kV. The smooth curve 
is a (normalized) fit to the data of the functional form 
(2vi/a 2 ) exp[— vf/a 2 ], with a = 54m/s, the stepped curve 
is the Monte-Carlo results for 5 kV. The negative values at 
high velocities are due to statistical noise in the data at short 
times. 



As a first measurement, the guiding signal was 
recorded as a function of the applied voltage. The re- 
sult, shown in Fig. ||(c), shows the expected quadratic 
dependence. To further characterize the source, the ris- 
ing slope of the signal, coming from the time-of-flight, 
was analyzed in more detail. By using fast switches, 
the turn-on time of the high voltage is less than 1 fis, 
much less than all other relevant time scales in the ex- 
periment. With the high voltage off, the density of slow 
molecules that could be guided decays rapidly with dis- 
tance from the nozzle. Therefore, it can be assumed that 
the molecules arriving at the detector must have entered 
the guide after the high voltage has been switched on. 
After a delay of a few ms, which depends on the applied 
voltage, the fastest guided molecules arrive at the QMS 



and the signal starts rising until it levels off after w 10 ms. 
From the delay and the rising slope, the longitudinal ve- 
locity distribution can be derived by differentiation. The 
analysis incorporates the following additional inputs : 1) 
the Monte Carlo simulation was used to estimate the 
velocity-dependent probability of hitting the detector af- 
ter emerging from the guide. This probability was used to 
perform a velocity-dependent correction, which accounts 
for the acceleration of the molecules when leaving the 
guide and the fact that longitudinally slow molecules are 
more likely to miss the detector, because they spread 
out over a larger solid angle; 2) making the reasonable 
assumption that the ionizing probability of the QMS de- 
pends on the molecular velocity and does not saturate at 
the lowest measured velocities, a velocity-dependent ion- 
ization probability was calculated. The resulting velocity 
distribution obtained for an electrode voltage of ±5 kV is 
shown in Fig. |j| Subtracting the signal with the shutter 
closed had only a small effect on the velocity distribution, 
and was omitted. 

Note that for a single molecular state the velocity dis- 
tribution should show a relatively sharp velocity cut-off. 
But given a mixture of states with different Stark shifts, 
the cut-off is smeared out. Surprisingly, this leads to 
a velocity distribution which can reasonably well be de- 
scribed by a 1-dimcnsional thermal distribution. For ex- 
ample, the smooth solid curve in Fig. represents a fit 
to the data with a temperature of 5.4 K. The experimen- 
tally determined velocity distribution is slightly narrower 
than the one from the Monte-Carlo simulation, and the 
mean velocity is shifted to a smaller value. This devi- 
ation is probably caused by experimental imperfections 
such as surface roughness of the electrodes, deviations 
from the design geometry, etc. The transverse velocity 
distribution of the guided molecules could not be mea- 
sured directly, however, simulations indicate a distribu- 
tion characterized by a temperature of 0.5 K inside the 
guide for ±5 kV electrode potential. 

The absolute flux was measured by calibrating the ion 
current from the QMS with the background signal caused 
by non-guided hot H2CO molecules. The partial pres- 
sure of the hot H2CO was measured with the QMS and 
related to an absolute pressure via a Penning pressure 
gauge [Balzers, compact full range gauge PKR260], tak- 
ing into account the relative ionization probabilities from 
Rcf. 1 14 . An ionization gauge [Varian UHV-24] con- 



firmed the Penning gauge readings to within 50%. The 
ionization volume is modelled by a (3 mm) 3 box, as in- 
dicated by the QMS supplier. Unfortunately, the exact 
dimensions of the ionization volume are hard to verify for 
our conditions. Therefore, we estimate a systematic er- 
ror in the absolute flux of the order of a factor of 2. The 
flux that was calibrated in this way is plotted in Fig. |] 
as a function of the absolute nozzle pressure. The mea- 
sured H2CO content in the nozzle is approximately 12%. 
In Fig. |4|, the nozzle pressure was deliberately increased 
to a point where the molecular flow approximation is no 
longer valid. As can be seen in the figure, for low nozzle 
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FIG. 4: The absolute measured flux through our guide as a 
function of nozzle pressure for an electrode voltage of ±5 kV. 
The symbols are data points, the curve a guide to the eye. 
For low pressures the gas in the nozzle is in the molecular flow 
regime and the flux is linear in the nozzle pressure. Above « 
0.02 mbar the mean free path for thermal molecules becomes 
less than the nozzle length and collisions start playing a role. 
As a consequence, the increase of the flux slows down and 
eventually the absolute guided flux is reduced. 



pressures the flux increases linearly, but starts deviating 
from the linear dependence above 0.05 mbar, reaches a 
maximum of nearly 10 9 s _1 , with peak density inside the 
guide of w 10 8 cm -3 , around 0.15 mbar before actually 
reducing with pressure. This behavior is caused by col- 
lisions. When increasing the pressure in the nozzle, the 
beam characteristics change from an effusive, thermal ve- 
locity distribution to a distribution which is peaked to 
a higher velocity. Moreover, slow molecules are kicked 
away by the fast ones in the guide and by background 
molecules. Therefore, it is important to keep the pres- 
sure in the nozzle sufficiently low. 



^From the Monte-Carlo simulation, the measured 
H2CO concentration and the theoretical flow impedance 
of the nozzle, the expected flux can be determined at 
a nozzle pressure of 0.025 mbar and an electrode volt- 
age of ±5kV. The result is 3.5 x 10 9 s~ 1 , a factor of 10 
larger than observed in the experiment. The discrepancy 
might be due to background collisions and, more likely, 
a non-perfect nozzle and a non-perfect alignment of the 
detector and the electrodes. 

To demonstrate the general character of the source, 
also a slow beam of deuterated ammonia (ND3) was pro- 
duced, with similar results to those obtained for H 2 CO. 
The deuterated species was used because of a null QMS 
background signal from contaminants at the mass peaks 
of interest and because its Stark shift is larger than that 
of NH 3 . 

In conclusion, we have experimentally demonstrated 
a continuous, high-flux source of translationally cold 
molecules at temperatures of a few kelvin. The molecules 
are conveniently delivered "on spot" , e.g. into an ultra- 
high vacuum, where they can be stored in an electrostatic 
trap or an electrostatic ring | ^5| . This simple and versa- 
tile method is generally applicable to all molecules with a 
reasonably large dipole moment. By reducing the radius 
of curvature of the bend, slower output beams can easily 
be achieved. In principle, the source can also be used 
as a novel kind of gas chromatograph that selects polar 
molecules from a non-polar buffer gas, or for enrichment 
of a gas with strong low-field-seeking states. 
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